GENERAL MIDPOINT SUBDIVISION 

QI CHEN AND HARTMUT PRAUTZSCH* 

Abstract. In this paper, we introduce two generalizations of midpoint subdivision and analyze 
the smoothness of the resulting subdivision surfaces at regular and extraordinary points. 

The smoothing operators used in midpoint and mid-edge subdivision connect the midpoints of 
adjacent faces or of adjacent edges, respectively. An arbitrary combination of these two operators 
and the refinement operator that splits each face with m vertices into m quadrilateral subfaces forms 
a general midpoint subdivision operator. We analyze the smoothness of the resulting subdivision 
surfaces by estimating the norm of a special second order difference scheme and by using established 
methods for analyzing midpoint subdivision. The surfaces are smooth at their regular points and 
they are also smooth at extraordinary points for a certain subclass of general midpoint subdivision 
schemes. 

Generalizing the smoothing rules of non general midpoint subdivision schemes around extraor- 
dinary and regular vertices or faces results in a class of subdivision schemes, which includes the 
Catmull-Clark algorithm with restricted parameters. We call these subdivision schemes generalized 
Catmull- Clark schemes and we analyze their smoothness properties. 

Key words, subdivision surfaces, midpoint subdivision, difference schemes, extraordinary 
points, characteristic map 
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1. Introduction. The Lanc-Ricscnfcld algorithm [LR80J for subdividing cardi- 
nal spline curves and tensor product surfaces has been fundamental to the devel- 
opment of box spline and general stationary subdivision schemes Pra84. CLR84, 
IDM841 IDGL911 ICDM91j . The algorithm is quite simple since it is based on itera- 
tively computing midpoints of the edges and faces of regular quadrilateral meshes. It 
can easily be generalized to arbitrary meshes and this is now called midpoint subdi- 
vision. It includes particular instances of the Catmull-Clark and Doo-Sabin schemes 
|CC78[ IDS78] , which even predate the Lane-Riesenfeld algorithm. 

A midpoint subdivision scheme consists of an operator A n ~ 1 R of degree n £ N, 
which is used successively to subdivide an input mesh M.. The refinement operator R 
maps M. to the quadrilateral mesh RM. by splitting every face of M. at its center as 



shown at the top of Figure 2.1 The averaging operator A maps M to the dual mesh 
AA4 that connects the centers of adjacent faces as shown in the middle of Figure [2~T] 
Lane and Riesenfcld developed their algorithm to generate known curves and sur- 
faces but the limiting surfaces of the midpoint subdivision schemes are not explicitly 
known, which complicates their analysis. In fact, Catmull had no proof for his scheme, 
which he presented in his dissertation |Cat74| and which he published later together 
with Clark |CC78j . The first analysis was outlined by Sabin in 1978 |DS78j but it 
took another 15 years until Reif established a rigorous proof |Rei93| . Reif's proof is 
based on an analysis of the subdivision matrix and requires (numerically) computing 
or estimating its spectral properties. The numerical nature of Reif's proof is a bot- 
tleneck when it is used to analyze infinite classes of subdivision schemes. Therefore, 
Zorin and Schroder were merely able to analyze the midpoint schemes up to degree 
9 |ZS01j and only recently further geometric arguments were developed that helped 
to show that all midpoint subdivision surfaces of any arbitrarily high degree are C 1 
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everywhere }PC11| . 

Many arguments used in |PC11| also apply to other subdivision schemes based 
on convex combinations and further generalizations of the results are possible. How- 
ever, certain technical difficulties and some fundamental challenges remain that are 
unsolved. In this paper, we address some of them and enlarge the class of subdivision 
schemes which can be analyzed, after appropriate extensions, by the technique devel- 
oped in IPCllj . Although we add to the averaging and refinement operators "only" 
the edge midpoint operator, which is simpler than the other two, we are challenged 
by new problems: 

• The limiting surfaces of regular meshes are only known for midpoint subdivi- 
sion and mid-edge subdivision. Again, if we had only finitely many schemes 
for which we could compute their symbols and run the standard analysis, this 
would not be a problem. In this paper, however, we consider infinitely many 
subdivision operators that can be decomposed into the three basic operators 
mentioned above. 

• The mid-edge operator changes the "orientation" of a mesh and therefore, we 
face additional topological difficulties. 

• The eigenvectors of a subdivision matrix have a basis of rotation symmetric 
eigenvectors of different frequencies. For midpoint subdivision, the eigenvec- 
tors with frequency are identical for all valencies of the extraordinary point 
but this is not the case for the general midpoint schemes we consider in this 
paper. This is particularly a problem in Section |I0| where we work with 
convex combinations other than midpoint averaging. 

General midpoint subdivision is defined in Section [2] and its C^-property for reg- 
ular meshes is shown in Sections [3] and [4] For arbitrary meshes, the subdivision is 
shown to be convergent in SectionjS] and its C^-property is analyzed in Sections [6j [71 
and[8j Some examples of general midpoint subdivision schemes are given in Section|9] 
and another generalization of midpoint subdivision, generalized Catmull- Clark subdi- 
vision, is defined and analyzed in Section [10] We provide a few concluding remarks 
in Section HT1 

2. General midpoint subdivision. First, we consider a generalization of mid- 
point subdivision including also the mid-edge or simplest subdivision operator V. This 
operator maps a mesh A4 to the mesh VM. that connects the midpoints of adjacent 
edges of A4 as illustrated at the bottom of Figure |2.I| 

More precisely, we consider the general midpoint subdivision operators 

U = A a >V v >R r °A a *- 1 V v °- 1 R r °- 1 ■ ■■A ai V Vl R ri 

for 

ai,ui,ri, . . . ,a s ,v s ,r s > 

with a + v > 1 and v + r > 1, where a ~ J2i=i a ii v — 12i=i v h an d r = 12i=i r *- 
The midpoint subdivision schemes A n ~ 1 R of degrees n > 2 and the mid-edge 

subdivision scheme V are special general midpoint subdivision schemes. 

We show that general midpoint subdivision surfaces are C 1 at their regular points 

and that they are also C 1 at their extraordinary points for a certain subclass of general 

midpoint subdivision schemes. 



QI CHEN AND HARTMUT PRAUTZSCH 



3 




3. A second order difference scheme. To analyze the smoothness of general 
midpoint subdivision schemes for regular (quadrilateral) meshes, we introduce a spe- 
cial second order difference scheme in this section and analyze its norm in the next 
section. 

A regular mesh C can be represented by the biinfinite matrix C = [cj]i G z2 of its 
vertices Ci, which are connected by the edges CjCj +ei , CjCj +e2 , j e Z 2 , where 



[d e 2 e 3 e 4 ] = 



1 1 
1 1 



To analyze smoothness, we need the (backward) differences 



1,2,3,4, 



and the mesh C, 



vv •" 



VVci 



of the second order differences 



VVCi := [ViViC; V 3 ViC; V 4 ViC; V 2 V 2 Ci V 3 V 2 Ci V 4 V 2 Ci] , 

where VjVfcCi = Vj(VfcCi). 

Let U be R, A, V, or a general midpoint subdivision operator. Since U maps a 



linear mesh [i(i)]iez 2 onto itself, where I : K 2 
order difference scheme Uy V exists such that 



is any linear function, a second 



(3.1) 



(UCU V = U^C 



vv"vv 
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for all meshes C (see |KobOO| Equations (11) and (12)]). 

Furthermore, for a mesh C = [ci] ieZ 2, each element of {VjVfcCj | j, k G {1, 2}, i £ 
Z 2 } is a linear combination of C^ v because 

ViV 2 C; = V 2 V lCi = V 3 V lCi - ViViC;^ . 

4. Smoothness for regular meshes. Applying a general midpoint subdivision 
operator 

U = A a °V v °R r ° ■■■A ai V Vl R ri 

to the grid Z 2 , we obtain the scaled grid Z 2 / 2 r+u / 2 or a rotated version of it, where 
r = r i ancl v — v i- We denote the scaling factor by 



(4.1) 



a = a{U) = 2- r ~ v ' 2 



and recall the following well-known fact, which can be derived easily from |Dyn92| 
Theorem 7.6]. 

Theorem 4.1 (Smoothness condition for regular control meshes). Let 

\\CWoa := sup ||ci|| and \\U\\ := sup \\UC\loo . 
iez 2 |c||oo=i 

If U maps 1? to crlr and if Us.— jo has some contractive power, i. e., if 



ui 

vv 



< a" 



for some k, then U is a C^-scheme, i. e., for any bounded mesh C there is a continu- 
ously differentiable function f(x,y) such that 



lim \\U k C~i (a k Z 2 )\ 



To check the prerequisites of this theorem, we use 
Lemma 4.2 (Estimates about the difference schemes). 



| A 



vvi 



1 



(a) 

(b) \\R^ V \\ =1/2, 

(c) \\V^\\ < 1/2, and 

(d) \\A^R^ V \\ <3/8. 

Proof. Since ^4y V = A and \\A\\ = 1, we obtain (a) 



Figure 4.1 shows a pair of meshes C and RC schematically. For these meshes, we 



have 



V 2 b 20 = V 2 b 21 = , 
V?b 30 = ^V 2 c 20 , 

V 2 b 31 = - (V 2 c 20 + V 2 c 21 ) , 

V 3 Vib 21 = ^V 2 ViCn = ~(V 3 Vic 21 - ViViC 21 ) , and 



V 3 Vxb 31 = - (V 3 V lC21 + V 2 c 20 ) 
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C = \c t 



BC = [b^-] 



Coo — boo bio c io — b2o b3o C20 — b4o 

Fig. 4.1. A pair of meshes C and RC. 




VC = [h i3 ] 



c 00 Cio c 2 o 

Fig. 4.2. A pair of meshes C and VC. 



Similarly, we can derive such equalities for all other elements in i?^ v Cy V . This proves 
(b). 

Using the notation shown schematically in Figure |4.2| for a pair of meshes C and 
VC, (c) follows from 

V|b 12 = V 3 V 2 b 12 = V 3 V 2 b 2 2 = ^V 3 Vic 2 i 
and from analogous equalities for all other elements in V^ V C^, V . 



C02 



-o 

Coi 



bio 

— — at h 



b i 
boo 



C12' 



1 

■ » ■ 



1 

b 22 1 ^ 

— 0- 

C21 



C22 



C = 



b2i 



1 

•t ■ 



1 



— 4. — 

1 



ARC = [by] 



C 00 T Tc 10 Tc 20 

Fig. 4.3. j4 pair of meshes C and ARC. 
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Finally, for the meshes C and ARC in Figure |4.3[ we have 
1 



V 2 b 21 



V 3 Vib 21 



16 
1 

16 
1 



(3V 2 c 21 



V 2 c 20 ) 



(3Vic 2 i + Vicu + Vic 20 - 5ViCi ) 



= (3V 3 V lC21 + V 4 ViCn + 2V 2 c 20 ) , 

VaVib,, = -^(V lC22 + Vicu + 3V lC21 - 3ViCu - 2ViCi ) 
16 

= i (V 3 V lC22 + V 4 V lCl2 + 2V 2 c 21 + 2V 3 V lC21 ) 
and analogous equalities for all other elements in A^, v i?^ v C^ v . Hence, 



l^vv^vvll ^ max 



3+1 3 



-1 + 2 1 



1 + 2 
16 



which establishes (d). □ 

Theorem 4.3 (C 1 continuity for regular meshes). Any general midpoint scheme 
U = A a 'V v 'R r " ■ ■■A ai V Vl R ri with a + v>landv + r>lisa C 1 -scheme for 
regular meshes, where a — a i> v = Sj= 



i i} and r = 

Proof. Using Equation (3.1) and Lemma 4.2 we obtain 



vvi 



4». 



< 

< 2 



2~2v—2r 

3 o-2r 

4 Z 
— u — It 



= a(U>) 



if w > 1 

if w = and a, r > 1 

2^ 



Since U 2 maps Z 2 to a (f/ 2 ) Z 2 , we conclude from Theorem 4.1 that C/ 2 is a C 1 - 
scheme. Since U and C/ 2 have the same limiting surfaces, U is also a C 1 -scheme. 
□ 

5. Basic observations. In this section, we consider general midpoint subdivi- 
sion for arbitrary (quadrilateral) meshes with extraordinary vertices or faces and we 
show that general midpoint subdivision converges, i.e., general midpoint subdivision 
surfaces are continuous. 

Interior faces and vertices of a (quadrilateral) mesh are called extraordinary if 
their valence does not equal 4. Subdividing by V, R, and A does not increase the 
number of extraordinary elements and isolates these elements. Therefore, it suffices to 



consider only (sub)meshes with one extraordinary element, as illustrated in Figure 5.1 
These meshes are called ringnets. 

Given a ringnet Af and a general midpoint subdivision operator 

U = A a *V v 'lC m A a - 1 V v - , -B?'- 1 ■ ■ ■ A ai V Vl R ri , 

we generate the sequence 

j^f(k) = jjkj^^ Furthermore, we subdivide just the regular 
parts of any N^* 1 and obtain for every k a limiting surface Sfc. Since is part of 
Sfc+i, it suffices to study the operator U 2 instead of U and we assume without loss of 
generality that v — v i i s even. Under this assumption, the "orientation" of UAf 

does not change from that of Af and the meshes AfW are either all primal (i. e., have 
no extraordinary face) or all dual (i. e., have no extraordinary vertex). 
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Fig. 5.1. Examples of rings and ringnets: a 1-ringnet with an extraordinary face of valence 5 
(left) and a 2-ringnet with an extraordinary vertex of valence 5 (right). The first rings Mi in both 
meshes are marked by bold lines. 



We say that a mesh Af influences another subdivided mesh M. if during subdi- 
vision every vertex in M has an effect on some vertex in M. and if additionally all 
vertices in M. depend on Af. 

Definition 5.1 (Ring, ringnet, core (mesh)). Let Af^ be the subnet of Af^ 
consisting of the extraordinary vertex or face of Af^ . The k-th ring around A/q is 
denoted by Aft- The mesh Af$...k consists of Nq and the next k rings of vertices 
around Afo and is called a fc-ringnet or fc-net for short. Furthermore, the submesh 
consists of Mi, . . . , Afj. 
The core or core mesh of Af (with respect to U) c onsi sts of all vertices influencing 



Afq for some k > 1. We denote it by Af c . Figure 5.2 shows some examples. The 
k-th ring around Af c is denoted by Af c .k- The subnet Af c .o...k consists of Af c and the 
next k rings of vertices around Af c and is called a c. /c-ringnet or c. fc-net for short. 
Furthermore, the submesh Af c .i...j consists of Af c .i, ■ ■ ■ , Af c .j if 1 < i < j ■ 



Fig. 5.2. The core meshes (in gray) and the rings Mc.l (bold edges) of a ringnet with an 
extraordinary face or vertex of valence 5 for the subdivision schemes VAV (left), AVAV (middle), 
and RV 2 R (right). The convex corners of Mc.l are marked by hollow dots. 

Depending on the context, we treat any mesh as a matrix whose rows represent 
the vertices or as the set of all vertices. 
It is straightforward to prove 

Lemma 5.2 (Dependence between nets after a subdivision step). 
(a) The subnet of AN , RAf , or VAf consisting of all vertices depending on any 
connected subnet of Af is connected. 
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(b) A/" c .o...fe determines A/" c g 2 k f or ^ — ^> e -> 

^.a.^fc = (C / AA c .o... fc ) c .o...2fe 



Lemma 5.3 (Dependence property of a core mesh). Let r > 0. Then, for 
some constant q depending on U and r, every vertex in Af c influences all vertices 
inM^ k l = ( uq+k A0c.o...r for all k>0, which is denoted by 



(g+k) 



Proof. Let a be such that A/"o... a contains Af c .o... r - For sufficiently large / and any 
k > 0, every vertex in A/" c influences all vertices in A/"o' +fe '', all vertices in 
. . . , and all vertices in + anc ^ hence also all vertices in M^ k ^ a \ Thus, we 

obtain the lemma with q = I + a. □ 

As explained above, every ringnet AA*) defines a surface s^. Since Sfe+i contains 
Sfe, we can consider the difference surface = Sfe_|_i\sfe whose control points are 
contained in a sufficiently large subnet A/^ q . with p not depending on k. Due to 



Lemma 5.2 (b), the operator U restricted to c.p- nets can be represented by a stochastic 
matrix S called the subdivision matrix, i.e., 



(5.1) 



KT (k+1) - SAf {k) 

- ,v c.0...p — ° - ,v c.0...p 



Theorem 5.4 (C -property of U). The subdivision surfaces generated by U are 
C° continuous. 

Proof. Since the subdivision matrix S is stochastic, i. e., S is a non-negative and 
real matrix and each row of S sums to 1, 1 is the dominant eigenvalue of S. Due to 



Lemma 5.3 there is an integer I > 1 such that 



This implies that S l has a positive column and, according to |MP89( Theorem 2.1], 
any sequence (S 1 c) converges to a multiple of the vector [1 . . . 1]' as i — > oo for all 
real vectors c. Therefore, the only dominant eigenvalue of S is 1 and it has algebraic 
multiplicity 1. 

Hence, the difference surfaces Sj\s,_i converge to a point and the surfaces gener- 
ated by U are continuous. □ 

To analyze the spectrum of the subdivision matrix S, we order any c.p-net J\f 
such that 

N c 

Mc2 



A/', 



c.p 



where A/" a consists of the convex corners and A4 of all other points in N c .i (see Fig- 



ure 5.2 for an illustration of the convex corners). With this arrangement, the subdi- 
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vision matrix S has the lower triangular form 



S = 



c 




* 


B 


* 


* A 


* 


* * 



where 

(5.2) 

(5.3) 

(5.4) 



A/;W = CN C , 



B 



and 



A 



To verify this, we recall that any point influencing the core mesh influences TVo and 
thus belongs to the core mesh. This implies Equation \b.2\ and shows that N c .\ influ- 
ences only points in {UM) c .i...oo- Hence, N c .2 influences only points in (UAf) c .2...oo, 
etc. We observe that Af a does not influence any point in Af^ and if U has a factor 



V, then Af a does not even influence Afa (see Figure 5.3). Thus, Equations (5.3 1 and 



(5.4) follow. Moreover, due to Lemma |5T2] (b), (UAf) c .2 is determined by (UAf) c .o...i, 
and {UN)c.3 is determined by (£W) c .o...2, etc. This implies the zero blocks in S. 

Hence, the eigenvalues of S are zero or are the eigenvalues of the blocks C, B, 
and A. 

Lemma 5.5 (Spectral radii of B and A). The spectral radii ps and pa of B and 
A satisfy 



PB < 2 



-r—a — v 



and pa < 



-r— a 



if v = 
if v > 



In particular, pb,Pa < °~{U) = 2 r v l 2 . 

Proof. Since B is non-negative, we get |HJ85t Corollary 6.1.5 on Page 346] 



PB < \\B\ 



1511 



where 1 



The vector Bl represents 7V C \ without the convex corners iiAf c — 0, Af — 1, Af a = 0, 



and N c .2...p = 0. Due to Lemma 5.2 (a), the vertices influenced by N c .i under A, R, 



and V form a connected subnet of AAf, RAT, and VAT respectively with an inner 
boundary whose vertices that do not form a convex corner have a value < 1/2. By 
induction, ps < 2~ r ~ a ~" follows. 

Similarly, we get p^ < A T T ~ a if v = 0. 

If v > 0, then decompose U into U{VUy. where U2 is the identity operator or a 



sequence of the operators A and R. Due to Lemma 5.2 (a), the vertices influenced by 
Af c .i form band rings with inner boundaries B\ in U2N and B2 in VU2M . For Af c = 0, 
we get that A/" c \ is determined by B2 which is determined by B\ excluding its convex 
corners (see Figure 5.3). Since the convex corners of N c .\ have no influence on B\ 
excluding its convex corners, we obtain pa = 0. □ 
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Bi 



B 2 



Fig. 5.3. B\ and B2 are the inner boundaries of the band rings in U%hf and in VU2M that are 
influenced by Afc.l- Note that the convex corner p of Bi has no influence on Bi- 



6. Symmetric ringnets. For the C 1 analysis of general midpoint subdivision, 
we need to investigate the eigenvectors and eigenvalues of the subdivision matrix 
S. We do this by subdividing special grid meshes as in [PCI lj and recall the basic 
definitions in this and the next section. 

Definition 6.1 (Grid mesh). A primal grid mesh of valence m and frequency f 
is a planar primal ringnet with the vertices 



K l '■■ = 



Mglj) 
Ms!,-) 



where g[ 3 = ie i2 * lf / m + j e ^(l+i)f/m £ c and i j>o i i e z m , i = 
A dual grid mesh of valence m and frequency f consists of the vertices 



hi 



1 



i.i-i 



g Z ■ 
01,3- 



+ sli), i,j>i,lez r , 



(see Figure 6.1). For fixed I, the vertices g' or with (i,j) ^ (0,0) of a grid mesh 



Af build the l-th segment ofAf. The segment angle of Af is Lp = 2itf/m. The half-line 
from the center gg through g' 10 is called the l-th spoke, denoted by Si(Af) or Si for 
short. 




Fig. 6.1. A primal grid mesh (left) and a dual grid mesh (right) with valence 5 and frequency 1. 



Topologically, any ringnet A4 is equivalent to a grid mesh Af. Therefore, we use 
the same indices for equivalent vertices and denote the vertices of A4 by p|,-. 

Definition 6.2 (Symmetric ringnet) . A planar ringnet of valence m with vertices 
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l 2 is called rotation symmetric with frequency f , if 

with 9 = 27r//m 



P' +1 



cos v — sin t 
sin 9 cos t 



A planar ringnet AT € K 2 is called reflection symmetric if its permutation Af consisting 



of the points p t 



(m-l)-l 



equals the conjugate ringnet Af consisting of the points 



J — 



r p l - ■ ' 




\ v l 1 

ri'i.X 






-/ 
L ^j,y J 



i. e., Af = Af . A rotation and reflection symmetric ringnet is called symmetric. 

7. The characteristic mesh. In this section, we construct a characteristic map 
of a general midpoint subdivision scheme 



U = A a -V v 'Fl!''A a - 1 V v -'-ET- 1 ■ ■ ■ A ai V Vl R r 



where Yfi 



Ei=i u i > !: E* 



■e: 



> 1, and Ei=i v i i s even - We follow 



PCll] and use results stated there for midpoint subdivision that are also valid for 
general midpoint subdivision since their proofs are only based on the properties that 
the subdivision scheme preserves symmetry and generates midpoints or any convex 
combinations. 

Theorem 7.1 (.Moo and A^). Let M. be the core mesh of a grid mesh with 
frequency f and segment angle 



if := 



2/tt 



€ (0,tt) 



whe 



denotes any matrix norm. Then the 



LetM k := (U k M) c / \\ (U k M) c 
following statements hold: 

(a) The sequence (A4k)keN converges to a symmetric eigennet M.^ with segment 
angle ip and a positive eigenvalue X v , which depends only on ip but not on 
f and m. (Moo)o...i has at least one non-zero edge. Additionally, we define 
Ajr = 1 7^1, where 7^ is the maximum eigenvalue associated with a rotation 
symmetric eigenvector with segment angle tt. 

(b) Restricting U to the core meshes, the eigenvalue X v is the dominant eigen- 
value of the eigenspaces of frequencies f and m — f and it has geometric and 
algebraic multiplicity 2. 

(c) X a > Xe > K for < a < 9 < tt. 

This can be proved as (5.4), (5.7), (6.3), and (6.4) in [PClT]. For midpoint 
subdivision schemes, is equal to 1/4 and to the subdominant eigenvalue [1q of 
frequency 0. This implies that A 27r / m is subdominant. However, for general midpoint 
subdivision schemes, A^ can be smaller than /j,q. Therefore, we use the following 
lemma to show that Xi^/m is subdominant for m > 5. 

Lemma 7.2 (Scaling factor a). 

(a) The scaling factor a = cr(U) defined in Equation (4-.1) is the subdominant 
eigenvalue of U for m = 4 and 

(b) A^ = a holds for any m and f such that = j. 

Proof. We consider a regular scalar- valued eigenmesh XAA = UA4 with eigenvalue 
A. Since 



X M 



vv 



vv M vv 
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and || (U 2 )^ v || < a (U 2 ) = a 2 (see the proof of Theorem 4.3 1, it follows that |A| < a 
or that = 0, meaning that M. is a linear image ot a regular grid Q, i.e., a 

linear combination of the constant mesh [1 ... 1]* with eigenvalue 1 and the two 
coordinates of Q . Since UQ = aQ, (a) follows for m = 4. Since there is a basis of 
rotation symmetric eigenmeshes, it suffices for m ^ 4 to consider a rotation symmetric 
mesh Ai with segment angle ir/2. Due to symmetry, the subdivided mesh UAA does 
not depend on /, whence (b) follows. □ 

This lemma together with Theorem |7.1| and Lemma |5.5| can be used as in the 
proof of Theorem (7.3) in |PCllj to derive the following corollary. 

Corollary 7.3 (Subdominant eigenvalue of U for valencies > 5). Let p be as 



in Equation (5.1), let U be a general midpoint subdivision operator mapping the space 
of c.p-ringnets of valence m to itself and let Ai be a c.p-grid mesh of valence m and 
frequency 1. Ifm>5, the meshes 

U k M 



\U k M\\ 



converge to a subdominant eigenmesh Moo of U called the characteristic mesh of U 
and its eigenvalue A 27r / m has geometric and algebraic multiplicity 2. 



Remark 7.4 (Subdominant eigenvalue for valence 3). Corollary 7.3 is also true 
for m — 3 if 

(7-1) A 27r / 3 > m&x{\p \,p B ,p A } , 

where ps and pa are the spectral radii defined in Lemma\5~. 



7.3 



Example 7.5 (The schemes VRV, VRVR, and V 2 for valence 3). Corollary 
holds for Ui = VRV and U2 — VRVR, because po = ps = Pa = 0, A 27r /3(C/i) = 1/8 
and \2tt/3(U2) — 1/16. However, for V 2 we get \2tt/3 = Po = Pb = 1/4 and hence 
C or pilar y\7 .3\ does not apply. 



8. Smoothness for irregular meshes. In this section we analyze general mid- 
point schemes of the form 

U = (A a °V v °R r °) ■ ■ ■ (A ai V Vl R ri ) =U 1 ---U n , 

where ^ a, + Vi > 1 < £ 2 v t + £i r, and U t € {A, R, V 2 , VAV, VA 2 V , . . .}. We 
call these schemes YASl-schemes and analyze their characteristic maps by estimating 
their partial derivatives using cones defined by direction vectors of the spokes So and 
Si of a grid mesh. 

To simplify the notation, we identify the real plane M 2 with the complex plane C 
by the bijection M 2 9 [x y]* 1 <-> x + ly <E C. 
Let 

C(a,0) := (o,oo)e J[min(a ' /5) ^ max(a ' /3)] 

and 

C o (a,0) :=C(a,f3)U{0} 

be the unpointed and pointed cones spanned by g 1111111 ^/ 9 ) an d e l max ( a ^) ] respec- 
tively. 

Moreover, we define an operator V 2 on a ringnet Af — [p^] by 
(8.1) V 2 (AA) := {V 2 p?, j - - p^-.x I * > Q,j > 0} . 
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Lemma 8.1 (A bound for V 2 ( C)). Let U be a YPN -scheme. Its characteristic 



mesh C = [cy] defined in Corollary 



7.3 and Remark 



7.4 satisfies 



m 2 



Proof. Let A4 be a sufficiently large grid mesh with valence to and frequency 1. 
First, we show that 



V 2 (U k M) C 2? := C 



2lT 7T 

TO ' 2 



holds for k > 0. Since V 2 (A^) C X>0) h suffies to prove that A, R, and VA l V map a 
symmetric ringnet A/" with V2(AT) C I> to a symmetric ringnet A/"' with V2(A/"') C XV 
Obviously, A/"' is symmetric. Furthermore, the elements of V2(AAT) and V2CRA/") 
are linear combinations of elements in V 2 (A/") with non- negative weights or, due to 
symmetry, are non-negative multiples of Ui or iuq, where Uo := 1 and Ui := e i i 2 ' n / m ) 
are the direction vectors of the spokes So and Si, respectively. For a net A l VAf, 
we consider the Ui -diagonals of all faces belonging partly or completely to the 0-th 
segment. The directions of these diagonals are either obtained by iteratively averaging 
the elements of VgJ^A /") or are, due to symmetry, parallel to Ui or m (see the top and 



middle of Figure 8.1 ). Hence, they lie in XV Since every element of V2(VA l VAf) is 
either obtained by halving the u 1 -dia gonals of A 1 VN or, due to symmetry, is parallel 



to zu (see the bottom of Figure 8.1 1, we conclude V 2 (VA l VJ\f) C X> . 

Second, we get V 2 (C) c X> because C = lim./.-^ U k M / \\U k M\\. 

Finally, we show V 2 (C) C V. Due to Theorem |7.1| (a), C0...1 has at least one 
non-zero edge. Let S be the subdivision matrix of U restricted to meshes of the same 



is 



size as C. For sufficiently large k, every element of V 2 (S k C) (= A fc V 2 (C), A > 0) 
a linear combination of V 2 (C) with non-negative weights and a positive weight for 
the non-zero edge in the 1-ringnet. Hence, V 2 (C) has no zero elements and thus, 
Va(C) CD. □ 

Theorem 8.2 (C^-property of VAV-schemes). Generic subdivision surfaces ob- 
tained by a VAV -scheme U are C 1 continuous at extraordinary points with valencies 



to > 5 and, if Inequality (7.1) is satisfied, also with valence m = 3 



Proof. Let c(x, y) : Q — ► C be three segments of the characteristic map of U, 



where = f2_i U Oo U as shown in Figure 8.2 and ciQ i is the i-th segment for 
i = -1,0,1. 

First, we observe that for a symmetric ringnet Af, each element of V 2 (AV), 
V 2 (2i?A/"), and V 2 (2yA'WV) is a convex combination of elements in V 2 (A/"), in 
V 2 (A") reflected at Si, and in — V 2 (A") reflected at So, where a reflected element 
has a weight which is less than or equal to that of the unreflected counterpart. Thus, 



by Lemma 8.1 and by induction, we see that V 2 (U k C/a k ) C T>q for k > 0. Since every 
partial derivative c y over f2o is the limit of a sequence of vectors Vfc £ V 2 (U k C/a k ), 
it follows that c y (Clo) C T> . 

Second, we show ^ c v (Qq). Any c !/ (x), x £ Oo, is a convex combination of £_i 
or £ , where Ei is the set of all y-edges in the segments i and i + 1 of U k C/a k for 
sufficiently large k. As in the proof of Lemma [8T| we observe 



£_i C C(2tt/to, tt - 2tt/to) D e 7r / 2 - 27r/m f , 
which implies ^ Cy(fio)- 



14 



GENERAL MIDPOINT SUBDIVISION 





Fig. 8.2. The domain f2 of c (left) and the —1, 0, 1-th segments of a grid mesh (right), where 
the y-edges in these three segments are marked by arrows and the y-edges in the 0-th segment are 
especially marked by double arrows. 



Hence, c y (£lo) C T> and similarly, c x (flo) C T> — tt/2. Since each element of T> is 



linearly independent with each clement of T>— 7r/2, c is regular over f2 (see Figure 8.3 1 
and hence, the total characteristic map of U is regular. Since c y (f2o) CD, c does not 
map any line segment between two points in to a closed curve, meaning that c is 



QI CHEN AND HARTMUT PRAUTZSCH 



15 



injective over fio- Moreover, since M. is a symmetric grid mesh whose 0-th segment 
lies in [0, oo) e^ 0, 27r / m l =: A and U preserves symmetry, it implies c(fi ) C A and c 
maps the interior of fio into the interior of A. Hence, the total characteristic map of 
U is injective. 




Fig. 8.3. The direction cones T> and T> — it/2. 



Thus, Reif's C 1 -criterion |Rei95[ Theorem 3.6] is satisfied, which finally concludes 
the theorem. □ 

Rem ark 8.3 (Checking the C 1 -property for valence 3). According to Theo- 
rem 8.2, a VAV -scheme generates C 1 surfaces at extraordinary points with valence 3 
if Inequality {7.1) is satisfied. Otherwise, if Inequality {7.1) is not satisfied, there is 
no guarantee that the eigenvalue A 27r /3 is a double subdominant eigenvalue and that 
it satisfies Reif's C -criterion as is the case with U — V 2 (see Example 7.5). In 
this case, we can use the more general C 1 -criterion in JPR081 Definition 5.10 and 
Theorem 5.25] with the weaker assumption that exactly two subdominant linearly in- 
dependent generalized eigenvectors exist while the characteristic map defined by the 
(non- generalized) subdominant eigenvectors is regular and injective. Regularity and 
injectivity still follow (even for valence 3 J if V 2(C) C T>o, where C is the control mesh 
of the characteristic map. 

9. Two examples. In this section, we consider two general midpoint subdivision 
schemes and study their smoothness properties. 

Example 9.1 (C 1 -property of V 2 ). Since V 2 is a YPN -scheme, it generates C 1 
subdivision surfaces for regular meshes and for meshes with extraordinary elements of 
valencies > 5 according to Theore ms\4-3\ and \8.2\ 

For valence 3, due to Example 7.5, Inequality {7.1) is not satisfied for V 2 . Accord- 
ing to \PR97f and \PR08[ Section 6. 3], the subdivision matrix of V 2 for valence 3 has 
exactly two linearly independent generalized eigenvectors associated with the subdom- 
inant eigenvalue 1 /4 and thus, the characteristic map exists. This map is regular and 
injective. This is shown in \PR97)j and can also easily be derived from Remark 8.3' 



The 0-th segment [cy] of the control net C = [cy] is given by 



c 13 
c° 

L 12 

C° 
c ll 



c 23 

C° 
L 22 

C° 
c 21 



L 33 
L 32 
L 31 



-1 1.75 

2.5 

1 3 



4 

4.75 

5 



V3i 



3 3.75 4 
2 2.5 2.75 
1 1 1 



Thus, the edges V2C°„- (J = 3,2,1, i = 1,2,3) of C are 



-l + V3i, -0.75 + 1.25 • -0.75 + 1.25 • y/St, 

-1 + s/3i, -0.5 + 1.5 • -0.25 + 1.75 • 

2 ■ 2 ■ \/3i, 2 • 



They all lie in T> — C(§, ^p) • Hence, due to Remark 
generates C 1 surfaces for valence 3. 



8.3, the scheme V 2 also 
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In general, if a general midpoint scheme is not a VAV-scheme, the C 1 analysis 
technique in this paper usually does not apply. 

Example 9.2 (Non VAV-schemes) . We consider the non YPN -scheme VRV R. 
It generates C 1 subdivision surfaces for regular meshes according to Theorem \^.'\ 

Operate this scheme on a dual grid mesh M. with segment angle ip. In Figure \9H\ 
the edge e in ( VRVR) M. has the direction angle tt — arctan(16 tan ip) and hence, it 
does not lie in the pointed cone T>q spanned by e lip and e l7r / 2 . 

Hence, for this scheme, we cannot check the regularity and the injectivity of the 
characteristic map by using the C 1 analysis technique in this paper (see the proof of 
Theorem 8.2). 




[Pi ' ' ' P4] 

[qi • • ■ qe] 
[ri • • • r 4 ] 

[si s 2 ] 



l+c 


5+c 


3+c 


l+3c 






2 


2 


2 


2 






s 
2 


s 
2 


s 
2 


3s 
2 


5 




2+c 
2 


5+2c 
4 


5+2c 
4 


3+c 
2 


6+3c 
4 


7+; 

4 


s 
4 





s 
2 


s 
4 


3 s 
4 


s 
2 



44 + 46c 24+8c 472+65c 24+8c 



32 
3s 
8 



16 

11.9 

32 



128 

3s 
8 



16 

13.9 

32 



676+257c 1352+515c 
512 1024 
51s 49s 
128 128 



Si - S 2 



-c/1024 s/64]* 



Fig. 9.1. Since V 2 (U 2 M) 3 e f V , V 2 (U k M/a k ) C V does not hold, where M is a dual 
grid mesh with segment angle ip, c = cos tp, s = sin tp, and U = VRVR . 



10. Generalized Catmull-Clark subdivision. In this section, we introduce 
the generalized Catmull-Clark subdivision schemes, which are another generalization 
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of the midpoint subdivision schemes, and analyze the smoothness of the subdivision 
surfaces. 

For an odd degree n = 2r + 1 > 3, we generalize the midpoint scheme 
M n = A n ~ 1 R = A 2 ■■■ A 2 R to M n =B r ■■■ B x R, 
and, for an even degree n = 2r + 2 > 2, we generalize the midpoint scheme 
M n = A n ~ 1 R = A A 2 ■■■ A 2 R to M n = AB r • • • B x R, 



where R is the refinement operator shown at the top of Figure |2.1| and where each 
smoothing operator Bi — B au p i is a generalized operator of A 2 with the mask shown 
in Figure 10.1 The functions cn{-) and ft(-) are non-negative functions depending on 
i and they satisfy < a* + ft < 1, a*(4) = 1/4, and ft (4) = 1/2. 



7/771 
/3/m 
7/m 



/3/m r)l/ m 




Bi — B au p. 

a := ai(m) 
P ■= ft(m) 
7 := 1 — a — /3 

a, [3 e [0,1) 

7 e (0,1) 



Fig. 10.1. Smoothing operator Bi 



g. for a vertex of valence m. 



It is straightforward to see that B x = A 2 holds on regular meshes and that 
M3 = B\ R is the Catmull-Clark scheme with restricted parameters a\ and ft, i. e., 
at extraordinary points with valencies 4), ot\ and ft satisfy ai(m), ft(m) € [0, 1) 
and ai(m) + ft(m) G (0, 1). 

If ai = 1/4 and ft = 1/2, then B\ = A 2 holds on arbitrary meshes and M n — 
B± R or M n = A B\ R is a midpoint subdivision scheme of degree n if n is odd or 
even, respectively. 

Theorem 10.1 ((^-property of M n ). LetM n =B r --- B X R orM n = B r ■ ■ ■ B X R, 
be a scheme as above with n > 2. Let A be the dominant eigenvalue associated with 
frequency 1 and let /io be the subdominant eigenvalue associated with frequency 0. If 
A > |/!o| for all valencies m > 5 or m = 3, then M n is a C 1 -scheme. 

Proof. First, for regular meshes, M n is the Lane-Riesenfeld scheme of degree n 
and thus, it is a C^-scheme. 

Seco nd, at extraordinary points, M n is convergent because, as in the proof of 



Theorem 



5.4 



it can easily be verified that the subdivision matrix of M n is stochastic 
with the simple dominant eigenvalue 1. 

Next, let M. be a rotation symmetric ringnet with non-zero frequency. If n is 
odd, we require M to be primal and otherwise to be dual. For any choice of the 
parameters on and ft in Bi, we get BiM = A 2 M for any primal rotation symmetric 
ringnet M with non-zero frequency and thus M n M. = M n M.. Hence, M n and M n 
share the same rotation symmetric eigenvectors with non-zero frequency and the same 
associated eigenvalues. Using Theorem 7.1 and the assumption X> \hq 
M n has the double subdominant eigenvalue A 



A 



Hence, M Tl 



we get that 
and M n have 



the same characteristic map. For the midpoint scheme M n , we have |/xo(Af„)| 
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A,r(M„) = 1/4 < A 27r / m (M„) and, according to Lemma 5.5 also /)s(M n ), pA{M n ) < 



1/4. Thus, Inequality (7.1) is satisfied and moreover, as in the proof of Theorem 8.2 
the characteristic map is regular and injective for all valencies m > 5 or m = 3. (One 
can also see [PClli (8.1) Theorem] for a proof of the regularity and injectivity of the 
characteristic map.) Thus, Reif's C 1 -criterion |Rei 95, Theorem 3.6] is satisfied, which 
concludes the theorem. □ 

Remark 10.2 (Generalized Catmull-Clark subdivision). The schemes M n can 
be further generalized by giving up the conditions a^(4) = 1/4 and (4) = 1/2. 
This generalized subdivision is called generalized Catmull-Clark subdivision and the 
resulting subdivision surfaces are analyzed in detail in the technical report [CP 12b], 
where we show that the generalized Catmull-Clark subdivision surfaces are C 1 at all 
regular points and at extraordinary points if one of the following conditions is satisfied: 

(1) m > 5 and cti, are constant functions. 

(2) m > 5, OLi,Pi are non- constant functions, and A > \no\. 



(3) m = 3 and Inequality (7.1) is satisfied, i. e., A 27r / m > max{|/i |, Pb-, Pa}- 
Here, m is the valence of the extraordinary element in the control mesh and A and po 
are defined in the above theorem. 

11. Conclusion. In this paper, two generalizations of midpoint subdivision, 
general midpoint subdivision and generalized Catmull-Clark subdivision, are intro- 
duced. They build two infinite classes of subdivision schemes, where the first class 
includes the midpoint subdivision schemes and the mid-edge, or simplest, subdivi- 
sion scheme and the second class includes the midpoint subdivision schemes and the 
Catmull-Clark subdivision scheme with restricted parameters. 

General midpoint subdivision surfaces for regular meshes are analyzed by esti- 
mating the norm of a special second order difference scheme and they are all C 1 
continuous. For irregular meshes, the smoothness is analyzed by a C 1 analysis tech- 
nique developed in |PC11| for midpoint subdivision. For this technique, we worked 
out several adaptions to the situation considered in this paper. In particular, the 
properties of the characteristic maps are analyzed without explicit knowledge of these 
characteristic maps and a deeper understanding of the spectral properties of the sub- 
division matrices at extraordinary points is provided. 

We see the results of this paper as a step towards even more general smoothness 
results in the spirit of the smoothness characterizations known for univariate corner 
cutting schemes |GQ96[ [PPS97J. In a consequent paper |CP12aj . the established C 1 
analysis technique is generalized to apply to specific infinite triangular subdivision 
schemes. We hope that our C 1 analysis technique can be applied to other subdivision 
schemes as well, e. g. schemes that can be factorized into general and simple convex 
combination operators. 
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